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Abstract 

The property of "return to equilibrium!^ is established for a class of 
quantum-mechanical models describing interactions of a (toy) atom 
with black-body radiation, or of a spin with a heat bath of scalar 
bosons, under the assumption that the interaction strength is suffi- 
ciently weak. For models describing the first class of systems, our 
upper bound on the interaction strength is independent of the tem- 
perature T, (with < T < To < do), while, for the spin-boson model, 
it tends to zero logarithmically, as T — * 0. Our result holds for inter- 
action form factors with physically realistic infrared behaviour. 

Three key ingredients of our analysis suitable concrete form 

of the Araki- Woods representation of the radiation field, Mourre's 
positive commutator method combined with a recent virial theorem, 
and a norm bound on the difference between the equilibrium states of 
the interacting and the non-interacting system (which, for the system 
of an atom coupled to black-body radiation, is valid for all tempera- 
tures T > 0, assuming only that the interaction strength is sufficiently 
weak) . 



*juerg@itp.phys.ethz.ch 

^present address: Department of Mathematics and Statistics, McGill University, 805 
Sherbrooke W., Montreal, Quebec, Canada, H3A 2K6, merkli@math.mcgill.ca 



1 Introduction 



The problem of return to equilibrium for models describing small systems 
with finitely many degrees of freedom coupled to a dispersive heat bath at 
positive temperature has been studied at various levels of mathematical preci- 
sion, since the early days of quantum theory. Fairly recently, a new approach 
to this problem based on spectral theory for thermal Hamiltonians, or Li- 
ouvillians, has been described and applied to simple models in [JP1, JP2]. 
The general strategy followed in our paper is based on the spectral approach 
proposed in these references; but our tactics are quite different and draw 
inspiration from techniques developed in [M] that have been motivated by 
methods in [BFSS]. For further results and methods relevant to our paper, 
see [BFS, DJ2, GGS1] and, in particular, [FM, FMS]. The work described 
in all these references relies on the deep insights of Haag, Hugenholtz and 
Winnink [HHW] and on the Araki- Woods representation [AW]. 

The main result proven in this paper is Theorem 2.1, asserting return to 
equilibrium for a class of models describing a "small system" with a finite- 
dimensional state space coupled to a "large system" , a dispersive heat bath 
at some temperature T, with < T < T < oo. The heat bath is modelled 
by a spatially infinitely extended free massless bosonic field. The systems 
we consider fall into two categories corresponding to a regular or a singular 
infrared behaviour of the coupling between the two subsystems. Both cases 
are physically realistic. 

We show return to equilibrium under the assumption that the interaction 
strength is sufficiently weak. For infrared-regular systems, such as toy atoms 
interacting with black-body radiation, our upper bound on the interaction 
strength only depends on T , but not on T < T . For infrared singular sys- 
tems, such as the usual spin-boson model, our upper bound on the interaction 
strength depends on T and tends to zero logarithmically, as T — > 0. 

The proof of Theorem 2.1, which is presented in Section 3, relies on a 
result of independent interest, Theorem 2.2, which says that the norm of the 
difference of the equilibrium states of the coupled system and of the non- 
interacting system is small, for weak interaction strength (uniformly in the 
temperature in the infrared-regular case, and with an explicit temperature- 
depent upper bound on the interaction strength for the singular case). The- 
orem 2.2 is proven in Section 4, and the proof draws on ideas developed in 
[A, FNV, BFS, F]. 
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With Theorems 2.1 and 2.2, we accomplish two goals. First, for infrared- 
regular systems, our results are uniform in the temperature T, for < T < T 
(where the high-temperature bound, T , has a clear physical origin, see also 
the comment after Theorem 2.1 below), assuming only that the interaction 
strength is small enough, with an upper bound only depending on T . Second, 
our results also hold for infrared- singular systems, provided the temperature 
is not too small (depending on the interaction strength). 

In order to render our discussion more concrete, we describe the mod- 
els studied in this paper more explicitly. The first class describes systems 
consisting of an atom, or of an array of finitely many atoms, coupled to the 
quantized electromagnetic field. We assume that the temperature T of the 
electromagnetic field is so small that it is justified to treat the atomic nucleus 
as static and to neglect the role played by atomic states of high energy, in 
particular those corresponding to the continuous energy spectrum describ- 
ing an ionized atom. Thus, the upper bound, T , on the temperature range 
considered in this paper is determined by the requirement that 



where k B is Boltzmann's constant, m at c 2 is the rest energy of an atom, and £ 
is the ionization (-threshold) energy. If T < T , with T satisfying (1), then an 
atom can be described, approximately (in the spirit of the Born-Oppenheimer 
approximation), as a quantum-mechanical system with a finite- dimensional 
state space spanned by those unperturbed atomic eigenstates corresponding 
to atomic energies E < const k B To in the discrete spectrum. This defines 
what we call a "toy (or truncated) atom" . 

The coupling of the toy atom to the quantized radiation field is described, 
in the dipole approximation, by an interaction term 



where e is the elementary electric charge, d at is the atomic dipole (moment) 

— * 

operator, and E is the quantized electric field. Furthermore, p is a density 
function corresponding to an approximate (smooth) 5-function peaked at 
the position of the nucleus and of width comparable to the size of the atom. 
(The interaction term (2) defines the Ritz Hamiltonian.) When expressed in 
terms of (Newton- Wigner) photon creation- and annihilation operators the 



k B T « m at c 2 , k B T < S, 



(1) 



— * — * 

ed at ■ E(p), 



(2) 
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interaction term (2) gives rise to a momentum-space form factor go(k) (see 
Section 1.1) corresponding to 



9o(k) = iy/\k\p{k) oc y/\k\, (3) 

for \k\ — > 0, where k is the photon momentum. Interactions characterized 
by an infrared behaviour go(k) oc \k\ p , as \k\ — > 0, with p > —1/2, are 
called infrared-regular. Nowhere in our analysis will the helicity of photons 
play an interesting role. The helicity- (polarization-) index will therefore be 
suppressed in our notation, and we shall think of the heat bath as being 
described by a scalar field (instead of a transverse vector field). 

The second class of models deals with systems of a quantum mechanical 
spin S, with S • S — s(s + 1) (and usually s = 1/2) coupled to a heat bath 
described in terms of a quantized, real, massless scalar field (p. Before the 
spin is coupled to the heat bath it exhibits precession around an external field 
B pointing in the z-direction. Its dynamics is generated by a Hamiltonian 

— * 

H spin = eS z , with e oc \B\. (4) 

The interactions of the impurity spin with the heat bath give rise to spin-flip 
processes described by an interaction term e.g. of the form 

gS x (f(p), (5) 

where g is a coupling constant, and p is a density function as described above. 
The bound, To, on the temperature range considered is determined by our 
desire not to take orbital excitations of the particle (an electron, neutron or 
atom in a dispersive medium, such as an insulator) carrying the impurity 

— * 

spin S into account. 

When <p is expressed in terms of (Newton- Wigner) creation- and anni- 
hilation operators the interaction term (5) gives rise to a momentum-space 
form factor g , with 

g (k) = ^Loc -L, (6) 



for | A; | — > 0, where k is the momentum of a scalar boson in the heat bath. 
Interactions characterized by an infrared behaviour (6) are called infrared- 
singular. 
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The physical interest of the second model, the spin-boson model, is some- 
what limited. But it has often been used to illustrate the phenomena of 
interest to us in this paper. 

A general class of model systems reminiscent of the ones just described is 
introduced, in a formal mathematical way, in Section 1.1 below. In the fol- 
lowing, we attempt to clarify what we mean by "return to equilibrium" . Let 
C d be the state space of the "small system" (the toy atom or impurity spin), 
and let B(C d ) denote the algebra of matrices acting on C d . Let 2U denote 
the algebra of Weyl operators over a suitably chosen space of one-boson test 
functions describing the quantum-mechanical degrees of freedom of the heat 
bath. The Weyl operators, which are exponentials of field operators smeared 
out with test functions, are bounded operators, and the algebra 2U they gen- 
erate is a C*-algebra. The kinematics of the composed system consisting of 
the "small system" and the heat bath is described by the C*-algebra 

21 = B(C d ) <g> 2U, (7) 

and its dynamics, in the Heisenberg picture, is given by a one-parameter 
group {a t }, with teR denoting time, of *automorphisms of 21. Before the 
small system is coupled to the heat bath, a t = a tt o is given by 

<x t ,o = a? <g> a{ (8) 

where af(A) = e ltHat Ae~ ttHat , A G B(C d ), is the Heisenberg-picture dynam- 
ics of an isolated toy atom, H at is its Hamiltonian, and where a{ describes 
the Heisenberg-picture dynamics of the heat bath. We choose {a/} to be 
the *automorphism group of 2U describing the dynamics of free, relativistic, 
massless bosons, such as photons (but, as announced, we shall suppress ref- 
erence to their helicity in our notation). 

Let LO a p and uj^ be the equilibrium states of the small system isolated from 
the heat bath, and of the free heat bath, respectively, at inverse temperature 
(5 = (ksT)^ 1 . Let H denote the Hilbert space of state vectors of the com- 
posed system obtained from the algebra 21 in (7) and the equilibrium state, 
up fi , given by 

up fi = uf <g> u f p , (9) 

before the small system is coupled to the heat bath, by applying the GNS 
construction. Furthermore let Jl^o G Tt denote the cyclic vector in 7i cor- 
responding to the state u>p t o, and let 7r@ be the GNS representation of 21 on 
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Ti,. Since oup^ is time-translation invariant, in the sense that wpfi(at t o(A)) — 
up t o(A), for all A G 21 and all times t G R, there is a selfadjoint operator, L , 
called thermal Hamiltonian or Liouvillian, acting on Ti with the properties 

-Kp{a tfi {A))=e aL ^p{A)e~ aL °, (10) 

for all A G 21, and 

L fy,o = 0. (11) 

In order to describe interactions between the small system and the heat bath 
at inverse temperature f3, one replaces the (unperturbed) Liouvillian L$ by 
an (interacting) Liouvillian L\, which is a selfadjoint operator on Ti given by 

L X = L + XIp, (12) 

where Ip is an operator on Ti determined by a formal interaction Hamiltonian, 
such as those in (2) or (5). The interaction Ip has the property that the 
dynamics generated by L\ defines a *automorphism group {(Jt,x} of the von 
Neumann algebra VJtp C B(Ti) obtained by taking the weak closure of the 
algebra 71"^ (21). This means that, for every operator A G VJlp and arbitrary 
t G WL, the operator 

a tiX (A) := e itLx Ae- itLx (13) 

belongs again to DJlp. (For a representation-independent way of introducing 
interactions between the small system and the heat bath, see e.g. [FM].) 
Following ideas in [A], one can prove that, for a large class of interactions Ip, 
there exists a vector Qp t \ G Ti. with the property that the state 

ivp, x (A) := (n fitX , A{lp, x ) , AeMp (14) 

is an equilibrium state for the interacting system, in the sense that it satisfies 
the Kubo-Martin-Schwinger (KMS) condition for the interacting dynamics 
on the von Neumann algebra VJlp, described by a t>X ] (see [HHW], or [JP2, 
BFS, DJP], for an explanation of these notions). The property of return to 
equilibrium means that the equilibrium state on Tip given by cup^ is dynam- 
ically stable, in the sense of the following definition. 

Definition. The system described by the von Neumann algebra DJlp and 
the time-evolution a t>x on Tip (a so-called VF*-dynamical system) has the 
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property of return to equilibrium iff, for an arbitrary normal state u on Tip 
(i.e., a state on Tip given by a density matrix on H) and an arbitrary operator 
A E Tip, 

\im cu(a t>x (A))=ujp, x (A), (15) 

or (more modestly) 

1 /"* 

lim - / ds oj(a s , x (A)) = ojp, x (A), (16) 

t-+OC t Jq 

(return to equilibrium in the sense of ergodic averages). 

The convergence in (15) and (16) follows from the KMS condition for 
(jjfax and certain spectral properties of the interacting Liouvillian, L x ; see 
e.g. [JP2,BFS]. Because cop jX is invariant under the time evolution a t , x , the 
interaction XIp in (12) can be chosen s.t. 

L x Qp, x = 0, (17) 

i.e., zero is an eigenvalue of L x . If zero is a simple eigenvalue of L x then, 
as a fairly easy consequence of the KMS condition and the von Neumann 
ergodic theorem, property (16) holds, and if the spectrum of L x is absolutely 
continuous, except for a simple eigenvalue at zero, then (15) holds (this is 
again easily from the KMS condition and the RAGE theorem. Let us also 
mention that if the kernel of L x is simple then L x does not have any nonzero 
eigenvalues, see e.g. [JP3]). 

The purpose of this paper is to exhibit a class of physically interesting 
interactions with the property that, for all f3, with (A; B T ) _1 = (3q < f3 < oo, 
return to equilibrium in the sense of ergodic averages, (16), holds, provided 
the coupling constant A is small enough, 

< |A| < A , 

where, for infrared-regular interactions, Ao only depends on /3 , while, for 
infrared-singular interactions, A — > logarithmically, as (3 — > oo; see Theo- 
rem 2.1. This result relies, in part, on the following result: Given any e > 0, 
there exists a positive constant Ai(e) and a choice of the phases of the vectors 
Qp t x and fig ;0 suc h that 

\\Vp, x - fyj i0 || < e, (18) 
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for all A, with |A| < Ai(e); in the infrared-regular case, the constant Ai(e) 
only depends on e, but is independent of (5, and it decays to zero as f3 — > oo 
for infrared-singular systems; see Theorem 2.2. 

A proof of return to equilibrium in the stronger sense (15), and uniformly 
in the temperature < T < T < oo has been obtained already in [BFS] and 
in [DJ2], under the infrared conditions go(k) ~ \k\ p , (\k\ ~ 0) for some p > 
and p > 2, respectively. In addition, [DJ2] show (15) in the infrared-singular 
case (6), for small coupling, tending to zero as T — > 0. The infrared condi- 
tions we impose to show (16) are p — —1/2 (T-dependent smallness of the 
coupling), and p — 1/2, 3/2, p > 2 (small coupling, uniformly in T). 

Acknowledgements. We thank V. Bach and I.M. Sigal for countless dis- 
cussions on related problems and spectral methods without which this work 
would never have been done. M.M. is grateful to V. Jaksic for illuminating 
discussions. We have enjoyed the hospitality of IHES during initial and final 
stages of this work. 



1.1 The model 

We consider a quantum system composed of a "small" subsystem interacting 
with a "large" subsystem. The pure states of the small subsystem, which is 
also called atom (or spin), are given by rays in the finite dimensional Hilbert 

space 

H at = C d . (19) 
The atomic Hamiltonian H at has simple eigenvalues E < E x < • • ■ < E d _i, 

H at = dmg(E ,E 1 ,...,E d _ 1 ). (20) 

It determines the dynamics af of observables A e B(H at ) according to 

af(A) = e itHat Ae~ itHat , (21) 

where tel. For any inverse temperature < (3 < oo there is a unique f3- 
KMS state on B(7i at ) associated with the dynamics (21), called the atomic 
Gibbs state (at inverse temperature f3). It is given by 

tr ( e - f3Hat ■) 

<(•) = L-W: • (22) 
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where the trace is taken over 7i a t- 

The large subsystem is infinitely extended and is described by a free, 
scalar, massless Bose field. Its state is taken to be the equilibrium state at 
inverse temperature < (3 < oo. The description of this state and the GNS 
representation is standard (see e.g. [AW], [JP1,2], [FM]). We present only 
the essentials and point out a modification we introduce (namely the phase 
in (36)). Let 

L\ := L 2 (R 3 , d 3 k) n L 2 (R 3 , W^k) (23) 

and denote by 2U(£q) the Weyl algebra over L 2 ,, i.e., the C*-algebra generated 
by Weyl operators W(f), f G L 2 ,, satisfying the CCR 

W(f)W(g) = e -$ Im{f ' 9) W{f + g) = e' llm ^W{g)W{f), (24) 

and the relations W(f)* = W(-f), W(0) = t (unitarity). The brackets 
(•, •) in (24) denote the inner product of L 2 (R 3 , d 3 k). The large subsystem is 
described by the (3-KMS state uj^ on 2U(Lq) associated with the dynamics 

a{(W(f)) = W(e u "f), (25) 

with dispersion relation 

u(k) = \k\. (26) 

An interaction between the two subsystems can be specified in a repre- 
sentation independent way in terms of a suitable *automorphism group a tt \ 
on the C*-algebra B(H a t) ®2U(-^o)> where A is a perturbation parameter and 
a t fl = af ® a{ . Here we do not discuss this procedure of defining a t> \ - this 
has been discussed in [FM]. Rather, we directly specify how the interacting 
dynamics acts (is implemented) on the GNS Hilbert space corresponding to 

cc^o = uf ® ujj, (27) 

the (/3, « t)0 )-KMS state on the algebra 21 = B{H at ) ® %0(L 2 ). The GNS 
representation of the algebra 21 determined by the state (27) is explicitly 
given in [AW] and has been put, in [JP1,2], in a form adapted to the use 
of the theory of spectral deformations (and of positive commutators). We 
use a slight modification of the representation in [JP1,2]. The representation 
Hilbert space is 

n = n at ® n at ® t, (28) 
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where 



T = T{L 2 (^L x S 2 ,du x da)) 



(29) 



is the bosonic Fock space over L 2 {R x S 2 ,du x da), where da denotes the 
uniform measure on S 2 . We use the following notational convention: we 
write L 2 (R x S 2 ) for L 2 (R x S 2 , du x da) and L 2 (R 3 ) stands for L 2 (M 3 , d 3 k), 
or for L 2 (IR + x S 2 ,u 2 du x rfcr) (polar coordinates). 
The cyclic vector representing in Ti is 



fy,o = off ® n. 

Here f2 is the vacuum vector in T and 



(30) 



QClt 



(tre-^r^e-^/V;®^, 



(31) 



where </?j is the eigenvector of iJ at associated to the eigenvalue Ej, see also 
(20). To complete our description of the GNS representation of (27) we need 
to give the representation map np : B(H a t) <8> 2^(-^o) — * £>(7~0- It is the 
product 

(32) 



with 



Tip = 7T 



7T at (A) 



.at 



(33) 
(34) 



and where, for h £ L 2 ( 
by 



x S* 2 ), </?(/i) is the selfadjoint operator on JF given 



<p(h) = 



a*(h) + a(h) 

71 



(35) 



The operators a*(h) and a(h) are standard creation and annihilation oper- 
ators on T, smeared out with h. We take h i— > a*(/i) to be linear. The 
real-linear map rg : L 2 , — > L 2 (M x S* 2 ) appearing in (34) acts as 



(rpf)(u,a) = ^ 



u 



u f(u, a), u > 0, 

e^J(-u,a), u<0, 



(36) 



where we represent / in polar coordinates and / means the complex conjugate 
of /. We have introduced an arbitrary phase <fi G M which can be chosen 
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conveniently so as to tune discontinuity properties of the r.h.s. in (36) at u — 
0. The origin of this freedom can be explained as follows. The expectation 
functional of is given by 



Ll 3 f 4{W{f)) = exp [-1 £ (l + -^-^ |iW| 2 ^ 



(37) 



which corresponds to the state of black body radiation at inverse temperature 
(3, see [AW]. We define a family of (equivalent) representations of the Weyl 
algebra 2B(£q) on the Hilbert space (29) by the map 



7T 



U+,U- 



(W(f)) 



exp 



lip T„ 



(38) 



where y2 is defined in (35), U+, U- are arbitrary unitary operators on L 2 (M 3 ), 
and 



It is easily seen that, for any choice of the unitaries U±, 



u(U + {l-e-^)-y*f_)(u,a), 
u(U-(e-P u - I)' 1 ' 2 f)(-u, a) 



u > 0, 
u < 0. 



(39) 



Q, exp 



lip 



n 



equals the r.h.s. of (37). Expression (39) reduces to (36) for 



id, 



U- = e i(t> . 



Remark. We recall that there is a second representation, 7ip +,U of QB(Lq) 
on T given by 



7T 



U+,U- 



(W(f)) 



exp 



tip 



■ u+,u - (e~ Pu/2 f 



(40) 



which commutes with the representation t^ + ' U ~ . 

In previous articles involving this setting, [JP1,2], [BFS], [DJ1], [M], 
[DJ2], [FM], [FMS], the freedom of choosing U± arbitrarily was not used, 
only U± — ±id was considered. For a suitable choice of U± one can apply 
the existing positive commutator methods, based on the generator of trans- 
lations in u G M as conjugate operator, to models with fermionic or bosonic 
fields having dispersion relation different from (26). These matters will be 
pursued in another work. Here we restrict our attention to the representation 
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(36), where is a phase determined by the interaction, see assumption (Al) 
and the discussion thereafter. 

We are now ready to define the interacting dynamics as the *automorphism 
group 

<TtA-) = e aL »(-)e-* L * (41) 
on the von Neumann algebra 

m p := 7cp (B(H at ) ® fflJ(Ljj))" C B(H), (42) 

where " denotes the double commutant (weak closure), and where the gen- 
erator L\, called the standard Liouvillian of the system, is the selfadjoint 
operator on H given by ([JP1,2], [FM]) 

L X = L + XI, (43) 

with 

L = L at + L f , L at = H at <g> l at - t at <g> H at , L f = dT(u). (44) 

Here, dr(-u) denotes the second quantization (acting on J 7 ) of the operator 
of multiplication by u G R, A is a coupling constant, and / is the finite sum 

/ = { Ga lat <PM9«)) - lat ® CatG^Ct ® ^(^(e"^ 72 ^))}, (45) 

where the operators G a are bounded, selfadjoint operators on Ti, at , and the 
functions g a G L\ are called /orm factors. C at is the antilinear operator of 
component-wise complex conjugation in the basis {<Pj}jZo diagonalizing H at . 
Note that L does not depend on the choice of the phases 0±, but I does. 
The following relative bounds are standard 

\\I(N+l)-^\\, \\(N+l)-^I\\<C(l + l/P), (46) 

where C is some constant which is independent of f3. 

At temperature zero {(5 = oo), the Liouvillian (43) corresponds to the 
Hamiltonian 

H x = H at + dT(iu) + \J2G a ® <p(g a ), (47) 

a 
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which describes interactions of the atom with the quantized field involving 
emission and absorption of field quanta. 

The pair (971^, a t ,\) is called a W^*-dynamical system. For A = the state 
on Wlp determined by f2g )0 is a (/?, (j t]0 )-KMS state. It is well known ([A], 
[FNV] , [BFS], [DJP]) that the vector ' 

n AA := Zp\ e-^+^n^o G W, (48) 

where is a normalization factor, and 7^ is obtained from I by dropping 
the second term in the sum (45), defines a (/3, cr^j-KMS state on QHg. 
Before stating our results we make two assumptions on the interaction. 

(Al) The form factors are given by g a (u,a) = u v lj a (u,a), where p takes 
one of the values —1/2, 1/2,3/2 or p > 2, and the g a satisfy a set of 
conditions we describe next. For fixed a and a, the map u \— > g a (u, a) 
is C 3 on (0, oo) and 



%9a\\i?(a*) < oo, for j = 0, 1, 2, 3. (49) 
If p = -1/2, 1/2 or 3/2 then the limits 

<&(0,a) := lim d J u g a (u,a) (50) 

exist, for j = 0, 1,2, and there is a phase O G not depending on 
a, a and j = 0, 1, 2, s.t. 

e-*&(0,a)£K. (51) 

In addition, if p = —1/2, 1/2 then we require d u g a (0,a) = 0. Finally, 
we assume that 

||« 2 0a||L3(R3) < oo- (52) 



(A2) It is assumed that 



min / da 

Em^En J g2 



^2 (<Pm,G a y n )g a (\E m - E n \,a) 



2 



> 0. (53) 
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Discussion of assumptions (Al) and (A2). Assumption (Al) concerns 
smoothness and decay properties of the form factors, which are necessary 
in the application of the Virial Theorem, see the remark after Theorem 3.2. 
If the interaction is characterized, according to (Al), by p = —1/2, then 
we choose the phase <fi in (36) to be = 20 o - For all other values of p we 
take = 7r + 20o- For p = —1/2, 1/2 an admissible infrared behaviour of 
the form factors is g a ~ u p times a constant, as u ~ 0. Other than for the 
applicability of the Virial Theorem, condition (52) is also used to show that 
L\ is selfadjoint (for any AgI). This follows from the Glimm-Jaffe- Nelson 
commutator theorem, see [FM]. 

Assumption (A2) is called the Fermi Golden Rule Condition and has been 
discussed extensively in previous works, see e.g. [JP1,2], [BFS], [DJ1], [M], 
[DJ2], [DJ3]. Its role is to guarantee that the probability of absorption and 
emission processes of field quanta does not vanish in second order perturba- 
tion theory (in A). This can be translated into a suitable positivity condition 
on an operator r , called the level shift operator, see (56) below. Let 



denote the projection onto the kernel of L , where P is the rank-rf projec- 
tion onto the kernel of L at , and Pq is the projection onto Cfl, Q being the 
vacuum vector in J 7 , see (29). We will see that if the non- negative operator 
UI5(L )IIl, where S is the Dirac distribution, has a one-dimensional kernel 
(the dimension is at least one, since the kernel contains the atomic Gibbs 
state (31)) then the system has the property of return to equilibrium. 

Theorem 1.1 Assume (53). There is an eo > 0, independent of f3 > (5q (for 
any (3q fixed), s.t. if < e < e then 



where C is a constant independent of (3, and Tq is a bounded operator on 
7i = Hat <8> 'Hat <8> J~, acting trivially on the last factor, J- ' , and leaving Ker L at 
invariant. Moreover, T Q restricted to Ker L at has zero as a simple eigenvalue, 
with the atomic Gibbs state as eigenvector, see (31), and is strictly posi- 



Tl = P ®P n 



(54) 



ui-^—m > r n - Ce 1 / 4 , 



(55) 
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Here, P Q at — 1 — P n at and P Q at is the projection onto Cf2^ . 

A proof of this result, in the case where the sum in (45) reduces to a single 
term, can be in found [M], [BFS]. It is easy to carry out that proof for the 
more general interaction (45). An explicit lower bound, in terms of (53), can 
be given: 

lip 

en 

which yields 70 in (56) by minimizing the r.h.s. over < (3 < 00. 



2 Main results 

Our main result on return to equilibrium is 

Theorem 2.1 (Return to equilibrium). Assume conditions (Al) and 
(A2). There is a constant Ao > ; independent of (5 > for any (5q > 
fixed, s.t. if 

r 1 ifp > -1/2 , 

0<|A|<A °{(l + lo g (l +i 3))^ ^=-1/2, (57) 

then the kernel of L\ is spanned by the interacting KMS vector Qp^x, (48). 
In other words, the system has the property of return to equilibrium. 

If the temperature of the heat bath is very large then second order processes 
of absorption and emission of field quanta do not dominate the ones of higher 
order, and we cannot expect to describe the physics of the system using per- 
turbation theory in second order (although, for the toy atom considered here, 
the phenomenon of return to equilibrium is expected to take place at all tem- 
peratures; see also our discussion in the introduction). This is why, in the 
following analysis, the dependence of the constant A in Theorem 2.1 on (3 is 
natural. The fact that, for p = —1/2, we must impose an upper bound on the 
coupling constant tending to zero, as (3 — > 00 (see (57)), can be understood 
as follows: Our methods are perturbative (in A) and rely on controlling the 
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(norm-) distance between the KMS states for the interacting and the non- 
interacting systems (see Theorem 2.2). One cannot, in general, expect this 
distance to be small, for small but non-zero coupling constants, uniformly 
in (3 — > oo. This is due to the fact that, for p = —1/2, and in the zero 
temperature limit, (3 — > oo, the groundstate of an interacting, infrared sin- 
gular system is not in Fock space (i.e., the Hamiltonian (47) does not have a 
groundstate in H at <g> ^(L 2 (M 3 )), see e.g. [AH], [Sp]), but the non-interacting 
system (A = 0) does have a groundstate in Fock space! Consequently, we 
expect the difference between the interacting and the non-interacting KMS 
state to diverge, as (3 — * oo, for p = — 1/2. 

Assuming that the interaction between the small system and the heat 
bath is such that (£lp,\, NQ^x) is small, for small values of A, uniformly in 
(3 — * oo, where N = dr(l) is the number operator in the Araki- Woods repre- 
sentation and Q/3,x is given in (48), then our methods can be used to establish 
return to equilibrium for sufficiently small values of A, uniformly in (3 > (3 , 
even when p = —1/2. 

From a more technical point of view, we can describe the above discus- 
sion as follows. A typical estimate involved in our analysis is inequality (46), 
where C is some finite constant. Then \\\I(N + 1) _1 ^ 2 || can be made smaller 
than any constant 5 > 0, provided |A| is chosen sufficiently small, indepen- 
dently of (3 > f3o, for an arbitrary, but fixed (3$ > 0. Similarly, in order to 
estimate the norm of the difference between the interacting and the non- 
interacting KMS state, we need an upper bound on the expectation value 
of the number operator N in the interacting KMS state Qp : x- As explained 
after the statement of Theorem 2.2, this expectation value is bounded above 
by \\\Ii(N + 1) _1 / 2 ||, where I\ is defined in (63). For p = —1/2, the latter 
norm is not uniformly bounded in (3 > (3q, but diverges logarithmically, as 
(3 — > oo. Thus, requiring it to be bounded by a small constant, we must 
assume that |A| \og{(3) is sufficiently small. 

Among the technical results used in our proof of Theorem 2.1 we single 
out the following one, which shows that the perturbed and unperturbed KMS 
states are close to each other, for small coupling constants. In the infrared- 
regular regime p > —1/2, the difference between the two KMS states is small 
independently of the inverse temperature. 

Theorem 2.2 Assume (Al) and let Pq^ x and Pn po denote the projections 
onto the spans of the interacting and non-interacting KMS states, Q/3 t \ (see 
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(48)) and Qp^ (see (30)), respectively. For any e > there is a A (e) > ; 
which does not depend on f3 > 0, s.t. if 

m , , x f 1 ifp > -1/2 . . 

HPn^-PnJI <e- (59) 

Remark. The constant A (e) in Theorem 2.2 depends on the spectral gap 
Ei — E > of the atomic Hamiltonian, and, if the norms are as- 

sumed to satisfy a ^-independent upper bound, then A (e) can be chosen 
independently of the dimension d of the atomic Hilbert space. 
We prove Theorem 2.2 in Section 4. 

3 Proof of Theorem 2.1 

We use a simplified version of the positive commutator (PC) method, intro- 
duced in the present context, for zero temperature systems, in [BFSS], and 
extended to the positive temperature situation in [M]. We refer to [DJ1], 
[DJ2], [O], [GGS1], [GGS2], and to the book [ABG], for recent different im- 
plementations of this method. 

3.1 Mechanism of the proof 

There are two key ingredients in our proof, the PC estimate and the Virial 
Theorem. While we give a proof of the PC estimate, we refer to [FM] for a 
proof of the Virial Theorem. 

Let Af = dT(id u ) be the second quantization of id u on T (c.f. (29)) and 

set 

A = i0\ (niR 2 t - R 2 JT^) , (60) 

where II = P <g> P n has been defined in (54), R e = HR € , TT = 11-11, 
R e = (Lq + e 2 ) -1 / 2 , and 9, e are positive parameters. We note that Aq is 
a bounded operator satisfying RanAo C T>(L\), and that the commutator 
[L, A Q ] extends to a bounded operator with 

ll[^,A ]||<C^ + ^. (61) 
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On the domain T>(N) of the number operator N = oT(]l) we define the 
operator 

B = N + \h+t[L x ,A ], (62) 

where 

h = J2 (G« ® it ® ^(3 M r^( 5a )) - i^C fl( G a C a ^^(9 u e- /3 ' l/2 r / 3(( ?(1 ))) . (63) 

a 

The operator B represents the quadratic form i[L\, Af + A ], see [FM]. 

Theorem 3.1 (Positive commutator estimate). Assume (Al) and (A2), 
and fix < r\ < 2/3. For any v > 1 set 

®„ = G ViN 1 ' 2 ) | HVH = 1, ||(AT + l) 1/2 Vll < 

There is a choice of the parameters e and 6, and a constant \i{rj) = Xi > 0, 
not depending on v and f3 > flo, s -t- if 

{ 1 ifp > -1/2 

° <|A|<Al minf 1 ^-" 2 \ if v- -1/2 (64) 

[ 111111 ^ l+log(l+/3) ' (l+log(l+/3))V lJP ~ 

then we have 

Pn , x BPn,, x > |A| 2 -V 3 - 9 "/ 2 7o P^ a , ( 65 ) 

in the sense of quadratic forms on R&nEA(L\) fl QS^, where A is any interval 
around the origin s.t. A fl a(L at ) = {0} ; E&{L\) is the spectral projection, 
and where 7 is given in (56). 

We note that it is enough, for our purposes, to examine B as a quadratic form 
on a subset of V(N 1 ^ 2 ), because any eigenvector tp\ of L x satisfies ip\ E 05^,,, 
for some u which is independent of |A| < 1. Moreover, for p > —1/2, u 
does not depend on (3 > (3 , while for p = — 1/2, u diverges logarithmically 
for large (3. These facts follow from the next result. 

Theorem 3.2 (Regularity of eigenvectors and Virial Theorem, [FM], 
[FMS]). Assume (Al). Letip\ be an eigenvector of L\. There is a constant 
c(p,/3) < oo, not depending on X, s.t. 

HA^aII <c(p,P)\X\U x \\, (66) 
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and s.t. for all (3 > fio (for any (3q > fixed), 

cM)<cM{ 1+ ^ (1 + p) filial (67) 

where C\ does not depend on (3 > f3 . Moreover, 

(B}^:=(^ x ,B^ x }=0. (68) 

Remarks. The constant c(p, (3) can be expressed in terms of the operator Ii 
given in (63) as follows: 

WhiN + iy^W <2^T||G Q || \\d u T P {g a )\\v = c(p,P). 

a 

One can understand (68) formally by expanding the commutator, 

tyx, [L\, A f + A }4> x ) = 2Tm (L x ip x , (A f + A )4> x ) = 0. (69) 
The same argument gives {[L X ,AA)^ = 0, from which it follows that 

> <A% - (A/iL x > (N)^-c(p,(3)\\\Ux\\\\N 1/2 i>x\\ 



fx 

1 1 



> ^(N)^--c(p,(3rX"\M\ 2 , (70) 

which yields the bound (66). In order to make the arguments leading to (68) 
rigorous, one needs to control multiple commutators of L x with Af + A of 
order up to three. In particular, we need the first, second and third com- 
mutator of / with the dilation generator Af to be a well-defined, relatively 
A rl / 2 -bounded operator, see [M], [FM]. The latter condition is satisfied pro- 
vided 

d J u Ti3(g a ) is continuous in u e M. for j — 0, 1, 2, and (71) 
dir p (g a ) G L 2 (R x S 2 ) for j = 0, 1, 2, 3. (72) 

We point out that for this argument, i.e. for the proof of (68), the L 2 -norms 
of the functions di u Tp(g a ) do not need to be bounded uniformly in (3. It is 
not difficult to verify that (71), (72) follow from (Al). Let p and <J)q De as 
in assumption (Al); then, for p — 1/2, 3/2, p > 2, we use the representation 
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(36) with = 20o, while for p = —1/2, we take = n + 20 o . 

The proof of Theorem 2.1 is an easy consequence of Theorems 3.1 and 
3.2. Indeed, if, for A satisfying (64), with v — v (introduced after Theorem 
2.1), there were a unit eigenvector ip\ e KerL A , orthogonal to f^A, then 

= (B)^ > |A| 2 -V 3 - 9 " /2 7o . (73) 

Relation (73) cannot hold since the r.h.s. is strictly positive. For p — —1/2 
condition (64) (with v = u = C[l + log(l + /?)]) gives (57), independently of 
V- 

3.2 Proof of Theorem 3.1 

Since £lp t \ is in the kernel of L\, the commutator B given in (62) cannot be 
strictly positive on the entire space; see (68). To show that dimKerLA = 1 
it is natural to try to show that 

B + 5P^ X > 7 , (74) 

for some 5 > 7, where 7 > 0. Let A C 1 be an interval around the 
origin not containing any non-zero eigenvalue of the atomic Liouvillian L at . 
In Subsection 3.2.1 we prove (74) in the sense of quadratic forms on the 
spectral subspace of L associated with the interval A (see (96)). Using this 
inequality, we show in Subsection 3.2.2 that 

P QptX BP a$t , > \lPn ^ (75) 

in the sense of quadratic forms on Rani^A'(-^A) H *B U , where E^(L\) is the 
spectral projection of L\ associated to an interval A', which can be chosen 
arbitrarily, as long as it is properly contained in A. 

3.2.1 PC estimate localized w.r.t. L 

We will use the Feshbach method with the decomposition 

H° A := Ran E° A = Ran E° A T1 © Ran ^TT, (76) 
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where II is given in (54), and where E A is the spectral projection of Lq 
associated with the interval A. For a presentation of this method resembling 
most closely the form in which it is used here we refer to [M], [FM], and, for 
more background, to [BFSS], [BFS], [DJ]. 

In what follows, C denotes a constant independent of A, 9, e, (3 > (3q (for 
any fixed /3q > 0), and C(p,/3) denotes a constant independent of X,9,e, 
satisfying the bound given in (67). The values of C, C{p,j3) can vary from 
expression to expression. 

From II = Po®Pn+Pn and the properties of A it follows that Ran -E^II C 
RanPn and 

E A Tl(B + 5P n ^)TlE A = E A TlN l/2 (1 + N' 1/2 \hN- l/2 ) N 1/2 UE A 

+E A U {i[L x ,A ]+5Pn, }X )UE A 

> ^E A U + E A m[L x ,A ]UE A 

> \(l-C 9 -p)EPjL, (77) 

provided 

\\P n N- 1 / 2 XI 1 N-^ 2 P n \\<C(p,P)\\\<l/2, (78) 
see the remark after Theorem 3.2, and where we use the bound 



Q\2 

\\E A n[L x ,A ]UE A \\<C— (79) 

which follows easily from the definition of A , (60). We choose the parameters 
s.t. 

C— < 1/2, (80) 



e 2 



and hence we have that 



E\ U(B + SPn 0A )UE° A > ^E° A U. (81) 



The Feshbach map associated with the decomposition (76) and with the 
spectral parameter m < 1/8, applied to the operator 

E A (B + 5P Q0 JE A (82) 
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viewed as an operator on the Hilbert space H A , is given by 

F n ,m{El{B + 6Pn jE° A ) = E° a u(b + 5P n ^ 

-(B + 6Pn gA )E° A n {B~+lPnZ-my 1 UE A (B + 5P n ^))uE A , (83) 

where the barred operator is understood to be restricted to the subspace 
RanZ^Il C H A . Using the definition of A , (60), and IUiiI = 0, one sees 
that 

npn = 29\ 2 niR 2 e m > o. (84) 

We show that the second term on the r.h.s. of (83), which is negative-definite, 
is smaller than UBU. By (81), the norm of the resolvent in (83) is bounded 
from above by 8 (for m < 1/8). Using this fact, the estimates ||L i2 e || < 1, 
Pell < e ~ 2 and ni[L A , A ]U = eXTiL x R 2 e m, we find that, for any ip e H° A , 
the modulus of the expectation value (-)^ = (ip, -tp) of the second term in 
the r.h.s. of (83) is bounded above by 

%\\E%n.{\h + i[L x , A ] + 5P^ X )U^\\ 2 

< i6^ 2 A 2 pjn^|| 2 

+c (V||np % , A n|| 2 + c( P ,p)\ 2 + 9 -^j \m 2 . (85) 

It follows that 

(F n , m (E A (B + 5Pn fj , x )E A ))^ 
> 26\ 2 (1 - 89) (UIR 2 £ IU^ + 5||P^ A n^>|| 2 

-c 6 ^ (lew) + ^ + ^Unp^nf) |H| 2 . (86) 

The expectation value on the r.h.s. of (86) is estimated from below using 

n/p 2 /n > i (r - c e 1/4 ) , (87) 

provided e < e , see (55), (56). Pick 9 and e s.t. 

9 < 1/16, e < e , (88) 
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and, for vp G Ran II, note the estimate 

6x2 Oft + 40 > — O^t + - Cel/A 



#A 2 



-7o 



(1 - ^ 1/4 /7o) Uf + ( g^Po,.* ~ p n 0fi/ ^ 



(89) 



where we use that P^att/; = Pn 0o ip for ip e Ran II. We choose 



0^2 ^2 

5 > — 7q > ^7o =: 7, (90) 



see also inequality (74), and 



e l/4 

C < 1/4. (91) 

7o 



The r.h.s. of (89) is bounded from below by 

^7o(3/4- ||P^ A -Pn„.oll)lM| 2 > ^7o U\\ 2 - (92) 

In the last step, we have applied Theorem 2.2, (59), which tells us that 
H-Pfy.A _ p %,oll < V 4 > provided 

A satisfies the condition (58) (with e = 1/4). (93) 
Combining this with (86), where we use 

Him^nii 2 = ||n(p n<M - P % , )nf < ||p % , A - p % , || 2 , 

gives 

a \ 2 

<Pn, m (P^(P + 5P % JP^)>^ > ^70 IHI 2 , (94) 

provided 

e„, „ x #A 2 e<5 2 



C^C(p,/3) + — + — j < 7o /4. (95) 
The isospectrality property of the Feshbach map tells us that 

El (B + 6P^ X )E A > min Q, ^70) E A = ^oE° A . (96) 
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3.2.2 PC estimate localized w.r.t. L\ 

Let < xa < 1 be a smooth function with support inside the interval A, 
s.t. Xa(0) = 1, and denote by x% — Xa(L ) and xa = Xa(L x ) the operators 
obtained from the spectral theorem. We show in this subsection that any 
unit vector ip e RanP^ x fl 53^, s.t. xa^> = i>, satisfies 

Q\2 

(B + SP^^iB)^— 7o , (97) 

provided suitable bounds on the parameters e, A, 6 are satisfied. We will 
repeatedly use the estimate 

II(i-x a)^II = II(xa-x^II < c|A| \\i(N + i)-^\\ \\(N + ly/^w 

< Cv\\\, (98) 

where the first inequality is a consequence of the standard functional calculus. 
Let us decompose the expectation value 

= (^(B + SPa^xl)^ (99) 
+ <(1 " X°a)(B + M> n<)i J(l - x°a)>^ (100) 
+2 Re <(1 - x°a)(B + SPa^)^ . (101) 

Because E° a x°a — Xa> inequality (96) implies that 

n\2 n\2 

(AiB + SP^JxDi > -^TollXA^II 2 > ^ 7 o(l-CV|A|)|H| 2 . (102) 
Since N + SPq p a is non-negative, we have that 

(100) > -|<(1-x a)(A/i+^^a,A ])(1-x°a)> 

> -|A| IKl-xi^ll \\h(N+l)-^\\ ||(7V + 1)V^|| 
-\\[Lx,A }\\ ||(1-Xa)^H 2 

n\2 / \2\ 

> -Cv 2 ^-[-C(p,P) + \\\ + ^y (103) 

where we have used (61). 

Our next task is to estimate (101). Since N commutes (strongly) with 
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Xa an d Pctp x ip = 0, and using that tp G RanP^ A , we conclude that 

Re<(l-XA)(5 + WV)X A }^ 
> S((l- Xa)P^ x (Xa - 1)L + ^ <(1 - Xa )(A7 1 + i[L A , A ]) XA ) t 



> -*||(1-Xa)^H 2 -C(p,)9)i/A s 



Xa)[Lx,A ] X ° a ), (104) 



Taking into account that (1— x A )n = and ||(1— % A )L 1 \\ < C (the constant 
is of the size | A| _1 ), one sees that the last term can be estimated as follows: 

((1-Xa)[Lx,A ] x ° a ) % 



= e\x\ 



l - x A )(A/n/i? £ - l x rju + \R e iui)x° A 



X 2 



ex 2 nx 



< Cv9\X\ — + |A| II^H 2 = Cv — + e 



(105) 



Plugging (105) into (104) and combining this with (102), (103), we arrive at 
the bound 



(B)+> 
9X 2 



4e 



Cv 



7o 



|A| 



—7o (1 - Ci/|A|) u-C(p, (3) + v\X\ + v— + ^ + e 



e 



(106) 

2 



Inequality (97) then follows by choosing parameters s.t. 



Cv 



7o 



A 2 



IAI 



Cv\X\ < 1/4 and — v»C(p, (3) + u\X\ + v— + ^ + e < 1/4. (107) 
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3.2.3 Choice of e, 9 and 5 

We must show that the conditions 

(78), (80), (88), (90), (91), (93), (95), (107) 
can be simultaneously satisfied. We set 
A = v- 9 / 2 X', 

some < e < 1, 



e = v- 6 \X 



3 1 \ / 1 e 



9 
6 



I*. 
9X 2 



some < t < e < 1 s.t. t > 3e - 2, 



-7o, 



(108) 

(109) 
(110) 
(111) 

(112) 
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and it is easily verified that there is a Ai > 0, depending on e,t, but not on 
v, (3 > Po, s.t. if 

< |A| < A x min (C(p, (3)' 1 , v 1/r ^ /2 C(p, (3)- 1/r >) , (113) 

where 77 = e — t > 0, then conditions (108) are met. The "gap of the positive 
commutator" (see (97)) is of size — = | A] 2_? V 3_9?? / 2 . The maximal value of 
rj under conditions (110), (111) is taken for e — > 2/3, t — > 0. 

4 Proof of Theorem 2.2 

The following high-temperature result is well known. Given any e > 0, there 
is an 77(e) > s.t. if 

(3\\\< 77(e) (114) 

then inequality (59) in Theorem 2.2 holds. A proof of this fact can be given 
by using the explicit expression (48) for the perturbed KMS state, and us- 
ing the Dyson series expansion to estimate \\£lp t \ — ^/3,o\\ (see e.g. [BFS]). 
Condition (114) comes from the fact that the term of order A™ in the Dyson 
series is given by an integral over an n-fold simplex of size (3, and, naively, 
(114) is needed to ensure that \\Q/3,\ — Q/3,o\\ is small. We shall improve our 
estimates on ||^,a — ^/3,o\\ by taking advantage of the decay in (imaginary) 
time of the field propagators. 

To start our analysis, we use the fact that the trace-norm majorizes the 
operator-norm to write 

ll-Pn^ - || 2 - H P ^/3,a - p n , o \\l = 2 (! - (tyj,A, -Pa^Aa)) 

< 2 (np tX ,p n atn^ + 2(n PtX ,p n Qp >x ), (115) 

where we use H — Pq < F^t + P n . Here, fijg* is the atomic Gibbs state at 
inverse temperature (3 given in (31), and Q is the vacuum vector in J 7 , see 
(29). We know that 

(n^Pn^x) < WN^n^W 2 < c(p,(3) 2 \X\ 2 , (116) 

where c(p, (3) satisfies (67), see Theorem 3.2. There is a (3\{e) > (3q s.t. if 
(3 > (3i(e) then 
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where ipo is the groundstate eigenvector of H at and Py ®y G B(H a t <8> 'Hat) 
is the projection onto the span of ipo <g> ipo. It follows from (115) that 

H^V - ^VH 2 < 2(^ x ,P Mo n l3 , x ) + e + 2c(p,/?) 2 |A| 2 , (118) 

for (3 > /3i(e). Let 

Q = G B(H at ) (119) 

be the projection onto the orthogonal complement of the groundstate sub- 
space of the atomic Hamiltonian H at so that 

< Q <g> l at + l at <g> Q. (120) 

Noticing that (fyg )A , Q <g> l at fi AA ) = (fi^, l at <g> Q = up >x (Q) we see 

from (118) that 

\\Pn,, x -PnJ\ 2 < 4^ x (Q)+e + 2c(p,(3) 2 \X\ 2 , (121) 

provided (3 > f3i(e). 

Proposition 4.1 For any e > i/iere exzsi ^(e) > and Ai(e) > swc/i 
£/ierf if P > fai^) and |A| < Ai(e) £/ien 

^,a(Q) < e- (122) 

The proof is presented below. For now, we use (122) to prove Theorem 2.2. 
We set 

(3 3 {e) := max(/3i(e), 02(e)), 

\' (e) := mm^X l (e),c(p,(3r 1 ^/7/2, V (e)/(3 3 (e)' ) ), (123) 

where 77(e) is the constant appearing in (114). In the case p > —1/2 the 
constant c(p, (3) has an upper bound which is uniform in f3 > (3q, see (67), 
and we take A (e) to be the r.h.s. of (123) with c(p, f3) replaced by this upper 
bound. For p = —1/2 we can find a A (e), indpendent of (3 > 0, satisfying 
(l + log(l+0))- 1 A o (e) < A'o(e), see (67). 

We always assume (58). Inequalities (121) and (122) yield 

\\Pn,, x - PnJ\ 2 < 6e, (124) 
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for > p 3 (e). If < p 3 (e) then (3\\\ < 77(e), and (59) follows from the 
high-temperature result mentioned above. This completes the proof of the 
theorem, given Proposition 4.1. 

Proof of Proposition It is convenient to work with a finite volume 

approximation 

of the KMS state up,\, where A = [-L/2, L/2] 3 C 1R 3 . (We introduce a finite 
box A just in order to be able to make use of some familiar inequalities for 
traces. The inequalities needed in our proof also hold in the thermodynamic 
limit, A /* M 3 ; but some readers may be less familiar with them.) In (125), 
the trace is taken over the Hilbert space H at <8> JF(L 2 (A, d 3 x)). For n = 
(ni,n 2 ,n 3 ) G Z 3 , let 

e*(x) = L- 3 / 2 e 2 ™*/ L , = 2 -^\n\ = + A + n 2 ) 1 / 2 (126) 



denote the eigenvectors and eigenvalues of the operator \J — A on L 2 (A, d 3 x) 
with periodic boundary conditions at <9A. We identify the basis {e^} of 
L 2 (A 3 ,d 3 x) with the canonical basis of / 2 (Z 3 ), and define the finite- volume 
Hamiltonian by 



if£ = H at + Hf + \v A , (127) 
5>«®¥>(fl£), (128) 



where c/ A e / 2 (Z 3 ) is given by 



«£<•») = (^J (r »^ 



and the operator 



L 



Hf = dT(hj), (130) 

acting on JF(/ 2 (Z 3 )), is the second quantization of the one-particle Hamilto- 
nian 

A A _ / EX, if (0,0,0), 
^ e «~i e A , if n = (0,0,0). 
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On the complement of the zero-mode subspace hf equals y— A with peri- 



odic boundary conditions. Changing the action of h 1 * on finitely many modes 



(always under the condition that e 1311 f is trace-class) does not affect the ther- 
modynamic limit. Similarly, we may modify the definition of g£ on finitely 
many modes without altering the thermodynamic limit. The existence of the 
thermodynamic limit, 

lim ujp X {A) — uj/3 t x(A), (132) 

can be proven by expanding e~^ H ^ into a Dyson (perturbation) series and 
using that 

<o(A) = l e _^ I (133) 

has the expected thermodynamic limit for quasi-local observables A. 

Our goal is to show that ujq X {Q) < e, for Q given in (119), provided (3 and 
A satisfy the conditions given in Proposition 4.1, uniformly in the size of A. 
In what follows, we will use the Holder and Peierls-Bogoliubov inequalities 
(see e.g. [S]). The Holder inequality (for traces) reads 

n 

pl..^! <nK-iu> ( i34 ) 

where 1 < pj < oo, J^jPj 1 = 1, and the norms are 

\\A\\ P = (tr \A\ p ) 1/p , for p < oo, and {{A^ = \\A\\ (operator norm). (135) 
The Peierls-Bogoliubov inequality says that 



tr (e A+B ) 



> exp [tr (Ae B ) /tre B ] , (136) 



tr e L 

which implies that 

< e P\H,^)\ = 1, (137) 



tr e~ pH o 
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since, by (128), uj^ {v a ) = 0. 

Using the Holder inequality one sees that, for any < r < (3/2, 

tr (e-^- 2 ^V^ A Qe-^ A ) 



< 



tr 



tre-^ A 



tr e 



tr <J [Qe'^^Q 



2M 



tr e" 



2M, for some MeN. 



-/JffA 



1 
2M 



, (138) 



where we are choosing r s.t. 

£- 

2r 
Setting 

= e - tH ov A e tH ° 

and using the Dyson series expansion we obtain 

Q e -^ H >Q = A + B, 
where the self adjoint operators A and B are given by 
A = Qe~^ H »Q 

B = E(-W 

n>l J ° 



(139) 
(140) 
(141) 
(142) 



Qv A (t n ) ■ ■ ■ v A (t 1 )e'^ H oQ dti • • • dt n . (143) 



'0<t„<...<tl< 7 



We plug (141) into (138), expand (A + B) 2M and use the Holder inequality 
to arrive at the bound 



<4AQ) < 



tr (\A\ 2M ) 

tre -/3//A 



+ 



tr (\B\ 2M ) 



(144) 



The first term on the right hand of (144) is easy to estimate. Let A = 
Ei — E > denote the spectral gap of the atomic Hamiltonian H at . Then 



tr (\A\™) tr nat (Qe-* g -*) 

tre~^< A ~" tr Wo( e-^»« 

<i-l 



d-l -P{Ej-E ) 



Ea—i 
J=l 6 



poc 

< V e^^-^) < 2 / e-^dx = 2 



/3A 



(145) 
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Taking into account (137) and (139), we obtain, for (5 > 1, 

"tr (\A\ 2M ) 



tre-^" 



i 

2M 



< 2e 



-2rA 



(146) 



In order to make the r.h.s. small, we take r large as compared to A -1 (hence 
(3 > 2r must be large enough). 

Next, we consider the second term on the r.h.s. of (144). From (137) one 
sees that 



tr (|-B| 2M ) a ( $Ha\ R|2aA _ , ,A /^/Wf£ r2m\ 



We expand 



where 



6\b^= £ TCfcx,...,^), 



(147) 
(148) 



fci,...,fc 2 M>l 



T(A; 1 ,...,A; 2M ) = (-A) fcl - 



2M 



dt\ 





(1) 



<*£> 



X 



2 _§_ +( 2 ) o j( 2M ) 

/ 2M d4 2) • • • / fc2_1 <> • • • / 4 2M) • • • / fc2M " 

•/o£r J(2M-1)JL J(2M-1) t 



dt 



(2M) 
&2M 



2M J ( 2M -!)W ^^Iffi 

xe^ A Q, A (0 • • V(ff>)Q g, A (4 2 2 } ) • • V(4 2) )Q x ■ ■ ■ 
..•xg/(C)-^f M) )^o. (149) 

Note that the time variables in the integrand are ordered, 

< t$ < ■ ■ ■ < t? < fg> < • • • < tf M) < (3. (150) 

Our goal is to obtain an upper bound on \up (T(ki, . . . , &2m))|, sharp enough 
to show that 

E Ko(m,---^2M))| (151) 

fcl,...,fe 2 M>l 

converges, and to estimate the value of the series. Note that the factors e@ H ° 
and e _/3H ° in the integrand in (149) drop when we apply (cyclicity of 
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the trace), and the expectation value of the integrand in the state Up = 
u f ® u p K ( see (133)) splits into a sum over products 

E ••• E ^(qg^.-.g^^q) 

(1) (1) (2M) (2M) \ 1 ' 

l ' ' ftj 1 ' ' k 2M 

^(^(©•••^ r ,(*i M °)), (152) 

where u/g* and u;^' A are the atomic and field KMS states at inverse tempera- 
ture (3, and 

G a (t) = e~ tHat G a e tHat (153) 
= e-^V(^)e^ A =a*(e-^^) + a(e^). (154) 

Using the Holder inequality (134) it is not difficult to see that 

2M 

QGatfV) ■ ■■G^ M) {tf M) )Q) < J] ||G>|| ■ ■ ■ \\G (155) 

Since u;^' A is a quasi-free state we can estimate the second factor in (152) 
with the help of Wick's theorem: 



, «(ti)---< JV (M)=E n , (ise) 

v (i,r)er 



where the sum extends over all contraction schemes, i.e., decompositions of 
{1, . . . , 2N} into N disjoint, ordered pairs (I, r), I < r. Applying (156) to 



4\<^ t{ S)---^ r At? M) )) (157) 



we find that all resulting terms can be organized in graphs Q, constructed in 
the following way. Partition the circle of circumference f3 into 2M segments 
(parametrized by the arc length) Aj = [(j — l)^g, J^m], j = 1, • • • , 2M. Put 

kj "dots" into the interval Aj, each dot representing a time variable t^ G Aj 
(increasing times are ordered according to increasing arc length). Pick any 
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dot in any interval and pair it with an arbitrary different dot in any interval. 
Then pick any unpaired dot (i.e., one not yet paired up) and pair it with any 
other unpaired dot. Continue this procedure until all dots in all intervals 
are paired; (notice that the total number of dots on the circle is even, as 
follows from the gauge-invariance of w^' A ) • The graph Q associated to such a 
pairing consists of all pairs - including multiplicity - of intervals (A, A') with 
the property that some dot in A is paired with some dot in A'. "Including 
multiplicity" means that if, say, three dots of A are paired with three dots 
in A', we understand that Q contains the pair (A, A') three times. The class 
of all pairings V leading to a given graph Q is denoted by Cg. Let 



(i,r)ev 



(158) 



denote the contribution to (156) corresponding to the pairing V. The nu- 
merical value, |£?|, corresponding to a graph Q is defined by 



Pec g 



and it follows from (156), (158) and (159) that 



(159) 



\ ki 1 



<£iei- 

Q 



(160) 



In order to give an upper bound on the r.h.s. of (160), we must estimate the 
imaginary-time propagators (two-point functions) 



g\ e -(P+ti-t r )h f 




A -(t r - tl )hj 



e 13 ^ - 1 



9i (161) 



e pn f - 1 

where the g A £ Z 2 (Z 3 ) are given in (129), and where U G A t , t r £ A r s.t. 
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< U < t r < (3. The r.h.s. of (161) equals 



n^(0,0,0) 



+g ai {2<iin/L)g ar {2<Kn/L)e-^- t ^ 



x 



e P^t - 1 



+ 



2tt 



[« 



+ e 



-(tr-tl) 



L J l ~ J - 1 ' 

In the limit L — > oo, the Riemann sum in (162) converges to 



(162) 



/ d 3 k IgaMga^e-^'-^ + g a ^k)g ar (k)e 



(t r -tt)\k\ 



oP\k\ 



J e 



p\k\ 



-, (163) 



since the form factors r satisfy conditions (Al). The term in (162) coming 
from n = (0, 0, 0) disappears in the limit L — > oo. (This shows why a 
redefinition of h A on the zero mode does not affect the thermodynamic limit.) 
It is not hard to see that, for arbitrary A;, A r and ti E A/, t r E A r , 



\U -t r \ > cL(A;,A r ) :-- 



P f o, 



if/ 



2M 1 \l - r\ - 1, lil^r 



and 



Defining 



/3- |t,-t r | > d+(A ; ,A r ) :=/3 



2M 



(|/-r| + l). 



d(A,A') := m in(d_(A,A / ),d + (A,A / )), 
we obtain from (161) and (163), and for L large enough, 



(164) 

(165) 
(166) 



f,A ( -t?H$ i A ■ 



e H11 fe r f 



<2 <C 



e -d(A,,A r )|fc| 



Given any two intervals A, A', set 



1/2 



-d(A h Ar)\k\ 
Va r i ^ _ e -/3\k\ y «r 



1/2 



(167) 



e -d(A,A')|fe| 

C(A, A') := 4 max ( g a , - fc| g c 



(168) 
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If L > C, for some constant C, then (168) is a volume-independent upper 
bound on the (finite-volume) two-point functions arising from contractions 
in the graph expansion (Wick theorem). We are now ready to give an upper 
bound on the r.h.s. of (160); (see also [F] for similar considerations). 

It is useful to start the procedure of pairing dots in the interval with the 
highest order k. Let n be a permutation of 2M objects, s.t. 

> K{2) >•••> K{2M)- (169) 

There are k.^i)) possibilities of pairing the dot t'f^ with some dot in an 

interval A (x(i)) . We associate to each such pairing the value 
'i 

fcjMi)) C(A 7r ( 1 ), A^ii))) < yj ^Jk^i)) C(A 7r (i), A^(x(i))), (170) 

where we use (169). Next, we pair the dot labelled by t^ (1)) (if it is still 
unpaired, otherwise we move to the next unpaired dot) with a dot in A (x(i)) 
and associate to this pairing the value 

\Jk-n(i) yjk^i)) C(A 7r (i), A^(i))). (171) 

We continue this procedure until all dots are paired. This yields the estimate 

2M 

Q j=i Q (A,A')eg 

Next, we establish an upper bound on the sum on the r.h.s. Using that 

\g a (k)\<C\k\^ (173) 

for some constant C, and for all a, provided \k\ is small enough, withp > — 1, 
it is easy to see that 

G(A ' Aj -°\ 1//3 + 1, d(A,A') = (174j 

Furthermore, using definition (166) and inequality (174), we see that, for any 
A, 

£ C(A, A') < T := C (l + i + -L. ( A) - 2 - 2 '\ < oo, (175) 
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provided p > — 1. Consequently, we find that 

l"j C(A, A') < T kl+ - +k2M . (176) 
g (A,A')ee 

Carrying out the integral over the simplex in (149), and using (152), (155), 
(160), (172), (176), we obtain the bound 



/ a \ki+-+k 2M 2M (u\kj/2 

\4 fl {T{k 1) ... ) k 2 M))\ < [c'w^) n^T-> ( m ) 

3=1 V 

13 \ki 1. 



where C = J2 a \\G a \\, and where the factor (^) 3 ^ is the volume of the 

simplex {t < t kj < ■ ■ ■ < ti < t + Thus, the series (151) converges for 

all values of A and (3 > 0, and 



Un B 2u\\ - c /| A | r Ay f r|A|rAV ( fc + 1 ) fc f , ( i 78 ) 

V /J - 1 1 2M^l 2M J (fc + 1) V ; 



Combining (144), (146), (147), (178), and using (139), we see that if L is 
large enough (independly of A or (3) then 

^ >A (Q) < 2e- 2 - A + C'lAlrr ^ (ClA|rr) fc . (179) 

The final step in the proof of Proposition 4.1 consists in showing that the r.h.s. 
(which is independent of A) can be made arbitrarily small, provided j3 is large 
enough and A is small enough. Pick (3 2 {e) > 1 so large that e~^ e ^ < e/2. 
For (3 > f3 2 (e) we choose r = (3 2 {e)/2 < (3/2. From the definition of T, (175), 
and the relation ^ = 2r, see (139), we see that Tr < C(e), uniformly in 
(3 > ^(e)- It follows that there is a Ai(e) > s.t. if |A| < Ai(e) then the 
second term on the r.h.s. of (179) is smaller than e/2. This completes the 
proof of Proposition 4.1. ■ 
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